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Abstract. The Gulf Coast of North America (GC) is a 'hot spot' of Phragmites diversity as several lineages (defined 
according to the haplotypes of their chloroplast DNA) differing in origin, genetic traits and phenotype co-exist and 
interbreed in this area. We analysed differences in photosynthetic characteristics among and within four haplotypes 
to understand if differences in gas exchange can be attributed to adaptations acquired in their native ranges. We 
collected rhizomes of four GC haplotypes (12, Ml, M and AI; including the phenotypes 'Land-type', 'Delta-type', 
'EU-type' and 'Greeny-type') and propagated them in a common controlled environment to compare photosyn- 
thesis-irradiance responses, C0 2 responses, chlorophyll fluorescence, the activity of ribulose-l,5-bisphosphate carb- 
oxylase/oxygenase (Rubisco), specific leaf area (SLA), pigment contents, stomatal density and guard cell length. The 
maximum light-saturated photosynthetic rate, A max , in the haplotype 12 (Land-type) and haplotype Ml (Delta-type) 
(34.3-36.1 jjimol C0 2 irT 2 s~ a ) was higher than that in the invasive Eurasian haplotype M (22.4 ± 2.3 ixmol 
C0 2 irT 2 s _1 ). The A max of haplotype AI (Greeny3-type) was 29.1 + 4.0 ixmol C0 2 rrT 2 s _1 and did not differ from 
the A max of the other haplotypes. The carboxylation rate (V cmax ) and electron transport rate (J max ) followed the 
same pattern as A max . The haplotypes also differed in SLA (17.0-24.3 m 2 kg -1 dry mass) and pigment content, 
whereas stomatal density and guard cell length, Rubisco activity and chlorophyll fluorescence did not differ signifi- 
cantly among haplotypes. The high photosynthetic activity and gas-exchange capacity of the two haplotypes origin- 
ating in tropical Africa and the Mediterranean area (haplotypes 12 and Ml) are apparently adaptations derived from 
their native ranges. Hence, the haplotypes can be regarded as ecotypes. However, it remains unclear how these dif- 
ferences relate to plant competitiveness and fitness in the GC of North America environment. 
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Introduction 

The common reed, Phragmites australis (Cav.) Trin. ex 
Steud, is a cosmopolitan emergent wetland grass occur- 
ring on all continents except Antarctica and in subtrop- 
ical to cold temperate climates (Brix 1999; Clevering 
and Lissner 1999). It grows well in wet habitats, along 
rivers and in the littoral zone of lakes down to ca. 3 m 
water depth due to its ability to aerate rhizomes 
and roots by efficient convective gas through-flow 
(Armstrong and Armstrong 1991; Brix et al. 1992; 
Armstrong et al. 1999). There is, however, considerable 
morphological variability between different populations 
and genotypes of P. australis, some of which may be 
explained by the large geographic range of the species 
and different ploidy levels (Clevering and Lissner 1999). 
The morphological, cytological and geographical vari- 
ation within the species accompanies substantial dif- 
ferences in ecophysiology and clonal development 
(Rolletschek et al. 1999; Kuhl and Zemlin 2000; 
Mozdzer and Zieman 2010; Achenbach et al. 2012; Eller 
and Brix 2012). Also, different genotypes of P. austra//s 
possess different degrees of phenotypic plasticity in 
photosynthesis when grown under common garden 
environments in different climates (Lessmann et al. 
2001; Eller and Brix 2012). 

In the Mississippi River Delta at the Gulf Coast of North 
America, five maternal lineages of Phragmites, defined 
according to their chloroplast DNA haplotypes (Saltonstall 
2002) and hereafter termed haplotypes, have been iden- 
tified (Hauber et al. 2011; Lambertini et al. 2012a). These 
haplotypes differ in their geographic origin and to a large 
extent in their phenotype. Nevertheless, there is no com- 
plete match between haplotype and phenotype, as one of 
the phenotypes belongs to three different haplotypes. 
This variation pattern is probably due to gene flow 
among haplotypes in the Gulf Coast, as evidenced by 
nuclear DNA variation (Lambertini et al. 2012a). 

The Land phenotype (or Land-type) is defined by its 
red, woody, branched stems and is prevalent inland 
from Texas to Florida, with only scattered occurrences 
in the Mississippi River Delta (Lambertini et al. 2012a). 
It is a hybrid between the species Phragmites mauritia- 
nus and P. australis, probably formed in tropical Africa 
where the two species co-occur. The Land-type has 
haplotype 12 like P. mauritianus in Uganda and Burkina 
Faso (Lambertini et al. 2012a). 

The Delta-type is the predominant phenotype in the 
delta and is identified by its tall shoots. It is an intro- 
duction from the populations of the Mediterranean 
region (Southern Europe, North Africa and the Middle 
East; Lambertini et al. 2012a) with which it shares 
haplotype Ml. 



The EU-type is identified by its purple inflorescence 
and small size compared with the Land- and Delta-types 
(also named 'short forms' by Hauber et al. 2011). It is a 
relatively recent arrival of the North American popula- 
tion of the invasive Eurasian P. australis with haplotype 
M which has spread in coastal wetlands along the Atlan- 
tic coast (Lambertini et al. 2012a). 

The Greeny-type occurs in small stands scattered in 
the Mississippi River Delta. The name refers to its charac- 
teristic blue-green stand colour. The Greeny phenotype 
belongs to three different haplotypes, namely haplotype 
M (hence Greenyl-type), haplotype AD (hence Greeny2- 
type) and haplotype Al (hence Greeny 3 -type). All Greeny 
haplotypes have been found in Europe (Lambertini et al. 
2012b), but one of them, the Greeny3-type, may be 
native to South Africa and introduced to the Mississippi 
River Delta from Europe. 

Despite the presence of haplotype M, the genetic 
pattern of the P. australis population at the Gulf Coast 
is very different from that of the East Coast (Lambertini 
et al. 2012a) where 14 closely related native American 
haplotypes and one distantly related Eurasian invasive 
haplotype (haplotype M) co-exist (Saltonstall 2002). In 
the Gulf Coast, three Eurasian haplotypes co-occur with 
one Mediterranean and one tropical African haplotype 
(Lambertini et al. 2012a). The European and Mediterra- 
nean haplotypes are closely related to each other and 
distantly related to the tropical African haplotype. In 
addition, haplotype M shows two distinct phenotypes 
in the Gulf Coast: the EU-type and the Greeny-type. 
While the EU-type is exclusive to haplotype M, the 
Greeny-type is shared by all European-related haplo- 
types (M, AD and Al). None of these haplotypes 
appears to be native to the Gulf Coast; haplotype 12 
with Land-type phenotype has been present in the Gulf 
Coast for the longest time (Lambertini et al. 2012a), 
whereas the introduction of haplotypes M, Ml, AD and 
Al probably occurred in the last couple of centuries 
(Chambers et al. 1999; Saltonstall 2002; Lambertini 
et al. 2012a). 

The aim of this study was to assess whether Phrag- 
mites haplotypes co-occurring in the Mississippi River 
Delta differ in their photosynthetic characteristics and 
whether these differences are shared by all genotypes 
of the same haplotype, hence being attributable to eco- 
types, or alternatively have resulted after colonization of 
the delta. We grew one to four genotypes of four haplo- 
types in a common controlled environment to ensure 
that possible differences in photosynthetic characteris- 
tics were due to genetic differences and were not 
affected by the variable environmental conditions in 
the delta. Achenbach et al. (2012) showed that ecophy- 
siological traits of P. australis are not related to 
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geographic range, but are genotype dependent. We 
therefore analysed intra-haplotypic variation in photo- 
synthetic characteristics in order to evaluate whether 
differences among haplotypes are consistent among 
genotypes within each haplotype or whether genotypes 
deviating in photosynthetic characteristics from the 
other genotypes of the same haplotype could be recog- 
nized by their phenotype (e.g. EU-type vs. Greenyl-type 
within haplotype M). 

Methods 

Plant material 

Eleven genotypes of P. australis were collected in June 
2009 from visually distinct stands in the Mississippi 
River Delta marshes (Table 1). Rhizomes from the geno- 
types were transferred to a greenhouse at Aarhus Uni- 
versity, Denmark (56°13'N; 10°07'E) and grown in 80-L 
tanks containing a mixture of sand and commercial 
compost. DNA was extracted from the genotypes for 
genetic analysis (Lambertini et al. 2012a). Haplotypes 
were previously determined by Lambertini et al. 
(2012a) based on chloroplast DNA sequences (tml-tml 
and rbcl-psal regions) by following the classification 
systems introduced by Saltonstall (2002). Phenotypes 
were determined based on nuclear DNA variation (micro- 
satellites, amplified fragment length polymorphisms and 
sequences), which matched the overall appearance of 
the stands in the field (Lambertini et al. 2012a). The geno- 
types were identified as 'haplotype 12, Land-type' (three 
genotypes), 'haplotype Ml, Delta-type' (four genotypes), 
'haplotype M, EU-type' (two genotypes), 'haplotype M, 
Greenyl-type' (one genotype) and 'haplotype Al, 
Greeny3-type' (one genotype). 

Rhizomes were transplanted to 6-L plastic pots con- 
taining a mixture of quartz sand and commercial 
compost, mostly consisting of Sphagnum (Pindstrup No. 
2, Pindstrup Mosebrug, Ryomgaard, Denmark). The pots 
were placed in a growth chamber (Bio 2000S, Weiss 
Umwelttechnik, Lindenstruth, Germany) under a day/ 
night cycle of 14/10 h, a temperature of 25/22 °C and a 
relative air humidity of 70/80 %. Light was provided by 
metal halide lamps at a photosynthetic photon flux 
density (PPFD) of ~150 |xmol rrT 2 s _1 at the base and 
500-600 ^mol m" 2 s 1 at the top of the plants. Each 
potted plant was placed in its own outer container to 
keep the soil water saturated. Plants were watered 
three times per week with a commercial fertilizer solu- 
tion (Pioner NPK Makro 10-4-25 + Mg and Pioner 
Mikro + Fe, Broste, Denmark). The genotypes were 
allowed to acclimate to the growth conditions in the 
growth chamber for 10 weeks before ecophysiological 
measurements were initiated. 



Photosynthetic light response 

The youngest fully developed leaves (the third or the 
fourth from the apex) of three shoots from each geno- 
type were used for the photosynthetic gas exchange 
measurements. Measurements were made with a 
Li-Cor 6400XT Portable Photosynthesis System (Li-Cor, 
Nebraska, USA) equipped with a Li-Cor 6400-02B LED 
light source. The sample chamber temperature was con- 
trolled at 25 °C and the relative air humidity at 35-50 % 
for all measurements. Photosynthetic light responses 
were determined from measurements at nine irrad- 
iances (2000, 1500, 1000, 700, 500, 250, 120, 60 and 
30 |jimol rrT 2 s~ a ) at 400 ppm C0 2 in the sample 
chamber. Leaves were acclimated in the leaf chamber 
for 3-5 min until steady-state gas exchange was 
achieved, and data then recorded using the light curve 
program of the Li-Cor IRGA. The light-response curve of 
each leaf was fitted using the quadratic equation of 
Prioul and Chartier (1977): 

A n = {<D,PPFD + A max - [(DiPPFD + A max ) 2 
-4/c$iPPFD] 05 )/(2/c)-R dark 

where A n (^mol C0 2 rrT 2 s~ a ) is the net assimilation 
rate; 3>j (mol C0 2 (mol photons) -1 ) is the initial slope of 
the light-response curve or apparent quantum yield; 
PPFD is the irradiance (|j,mol rrT 2 s~ a ); A max (|j,mol 
C0 2 rrT 2 s _1 ) is the light-saturated rate of gross photo- 
synthesis; k is the convexity; and R dar k (ixmol 
C0 2 rrT 2 s _1 ) is the dark respiration rate. The dark respir- 
ation (Rdark). the light compensation point (I c ), the light 
saturation point (I k ), and <t>, and A max were estimated 
from individual light-response curves using Photosyn 
Assistant software, version 1.1.2 (Dundee Scientific, 
Dundee, UK). The A n at 0 PPFD was inferred from the 
above equation. The data from the light-response 
curves at light saturation were used to calculate the in- 
trinsic water-use efficiency (WUE,) as the ratio of A max to 
stomatal conductance (g s ), also measured by the Li-Cor 
6400XT. 

Photosynthetic C0 2 response 

The leaves used for the light-response curves were 
marked and also used for C0 2 -response curves. The C0 2 - 
response curves were taken at nine C0 2 concentrations 
(800, 700, 600, 500, 400, 300, 200, 100 and 50 ppm). 
C0 2 was supplied with a 12-g C0 2 cartridge mounted 
in the Li-Cor 6400XT and assimilation was measured at 
a PPFD of 2000 |xmol itT 2 s 1 and at 25 °C. 

The carbon fixation kinetic models for terrestrial C3 
plants (Farquhar et al. 1980) as modified by Dubois 
et al. (2007) and Sharkey et al. (2007) were used to 
describe the relationship between the net carbon 
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Table 1. List of Phragmites genotypes used in the study. 



Genotype" 


Haplotype b 


GenBank accession no. 


Phenotype 0 


Distribution and origin 0 


ROMS3 


12 


HQ664450 + AY016334 


Land-type 


Distributed along the Gulf Coast of North America from Texas to 


ROM16 


12 


HQ664450 + AY016334 


Land-type 


California and scattered in the Mississippi Delta; originating in 


VVmjj 


1/ 


nyDDHHjU -f- HTUlDJJ't 


Land-type 


tropical Africa and anciently established in America 


R0MS4 


Ml 


JF271678 + AY016335 


Delta-type 


Dominant phenotype in the Mississippi River Delta; originating from 


R0MS7 d 


Ml 


JF271678 + AY016335 


Delta-type 


the Mediterranean region (North Africa, Middle East, Southern 


WHS4 


Ml 


JF271678 + AY016335 


Delta-type 


Europe) 


R0M4 


Ml 


JF271678 + AY016335 


Delta-type 




R0M2 


M 


AY016327 + AY016335 


EU-type 


Distributed throughout North America; distinguished from Land and 


WHS2 


M 


AY016327 + AY016335 


EU-type 


Delta phenotype by its purple panicles and smaller size; introduced 


0CT1 


M 


AY016327 + AY016335 


Greenyl-type 


from Eurasia 


Pal07gcUS 


AI 


AY016326 + HQ664451 


Greeny3-type 


Exclusively found in the Mississippi Delta; can be distinguished by its 
bright blue-green leaves, introduced from Europe 



"Labels according to Hauber et at. (2011). ROM, Romere Pass and its distributaries; ROMS, south of Romere Pass; WHS, White Splay; OCT, Octave Pass. 
b Haplotype 12 is the Gulf-Coast lineage (Saltonstall et at. 2004) and the cp-microsatellite variant 12 of haplotype I (Lambertini et at. 2012b). 
"According to Lambertini et at. (2012a). 

d Genotype ROMS7 (Delta-type) shares alleles with the EU-type genotypes. All genotypes are from the interior marshes of the Mississippi Delta. All the rhizomes 
were collected in June 2009. 



assimilation rate, A (jxmol C0 2 rrT 2 s~ a ), and the inter- 
cellular C0 2 concentration, O, (ixmol mor 1 ), and to esti- 
mate the key biochemical limitations to steady-state C3 
photosynthesis: 

A = min(A c ,/Aj, A p ) 
A c = [V cmax (q - r*)]/[q + K c (l + 0/K 0 )] - R d 
Ai = [JmaxCQ - r*)]/(4Ci + 8P) - R d 



275.7 mmol mor 1 ) and the photorespiratory C0 2 com- 
pensation point (r* = 28.7 ixmol mor 1 ) were adjusted 
using the temperature coefficients from Bernacchi 
et al. (2003) and Dubois et al. (2007). 

Carbon dioxide assimilation was assumed to be 
Rubisco limited at O, < 250 ixmol mor 1 ; hence, the A/O, 
datasets with O, < 250 ixmol mor 1 were fitted to A c to 
estimate V cmax and R d (Dubois et al. 2007). Then R d esti- 
mated from A c was used to estimate J max from A s using 
A/0, datasets with C, > 250 ixmol C0 2 mor 1 . Our data 
did not allow estimation of TPU. 



A p = 3TPU - R d 

In these equations, A c , Aj and A p are the net C0 2 assimi- 
lation rates (ixmol C0 2 nrT 2 s _1 ) limited by Rubisco, 
ribulose-l,5-bisphosphate (RuBP) regeneration and 
triose phosphate use (TPU), respectively. V cm ax is 
the maximum carboxylation rate of Rubisco 
(ixmol rrT 2 s~ a ); R d is the rate of mitochondrial respir- 
ation (ixmol C0 2 rrT 2 s _1 ); K c and K a are the Rubisco 
Michaelis-Menten constants for C0 2 (iximolmor 1 ) and 
0 2 (mmol mor 1 ), respectively; O is the partial pressure 
of oxygen (mmol mor 1 ); J max is the maximum rate of 
electron transport to reduce NADP+ for RuBP regener- 
ation (ixmol e~ nrT 2 s _1 ); and r* is the photorespiratory 
C0 2 compensation point (ixmol mor 1 ). The Michaelis- 
Menten constants (K c = 406.8 ixmol mor 1 , K a = 



Stomatal density and length 

Stomatal density (mm -2 ) and the length of guard cells 
(|xm) were determined on the third or fourth youngest 
leaf. Three new leaves produced during the 10-week ex- 
periment in the growth chamber were harvested from 
each genotype and fixed in 70 % ethanol for 7 days 
and then hydrated by placing them in 70, 50, 25 and 
0 % ethanol sequentially, for at least 15 min in each so- 
lution. The leaves were cleared in fuchsin-KOH for 24 h at 
60 °C, then rinsed in water and dehydrated by placing 
the leaves in water first and then sequentially in 70, 96 
and 99 % ethanol. Finally, leaves were embedded in 
Euparal (Carl Roth Gmbh, Karlsruhe, Germany). Stomatal 
density and the length of guard cells on both upper 
(adaxial) and lower (abaxial) sides of leaves were 
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determined under x40 magnification with a calibrated 
reticule in a light-transmission microscope. 

Specific leaf area and pigment analysis 

After taking photosynthetic light and C0 2 response mea- 
surements, the leaves were cut and their one-sided area 
was obtained by the weight/area ratio of photocopies of 
the leaves. The leaves were then lyophilized in a freeze- 
drier for the determination of dry mass (DM). The specific 
leaf area (SLA, m 2 kg -1 DM) was calculated as leaf area 
per dry mass. Pigments (mg g _1 DM) were extracted in 
96 % ethanol and the contents of chl a , chl b , chl a+b and 
total carotenes and xanthophylls (total carotenoids) 
were analysed spectrophotometrically according to 
Lichtenthaler (1987). Pigments (mg rrT 2 ) were obtained 
as the ratio between pigments (mg g~ a DM) and SLA 
(m 2 kg -1 DM) obtained as described above. 

Rubisco activity 

Three leaves of each genotype were harvested for the 
analysis of initial Rubisco activity (non-activated 
Rubisco) and total Rubisco activity following Hansen 
et al. (2007). Harvesting took place in the light and the 
entire leaf was immediately frozen in liquid nitrogen 
after excision. The frozen leaf was ground in a mortar 
containing liquid nitrogen. Approximately half a labora- 
tory spoon of the ground leaf material was then trans- 
ferred to a chilled mortar and further ground in 5 mL 
of extraction buffer containing 50 mM Bicine (pH 8), 
1 mM EDTA-Na 2 , 10 mM MgCI 2 , 5 mM dithiothreitol 
(DTT), 10 mM isoascorbate and 2 % (w/v) polyvinylpyr- 
rolidone. Initial Rubisco activity was determined in an 
assay solution consisting of 500 mM Bicine (pH 8), 
1 mM EDTA-Na 2 , 100 mM MgCl 2 , 50 mM DTT, 192.5 mM 
NaH 14 C0 3 and 5 mM RuBP. The reaction was initiated 
by adding ground plant extract and stopped after 60 s 
with 6 M HCl. Total Rubisco activity was analysed by 
activation in an assay solution for 5 min before 
addition of the 5 mM RuBP. The reaction was stopped 
after 60 s with 6 M HCl. Assays were carried out at 
25 °C in a total volume of 300 (il using 6-mL vials. 
Extracted samples from the assay were dried at 60 °C 
for 24 h and then re-dissolved in two drops of 6 M 
NaOH and 1.2 mL of ultra-filtered water. The amount 
of radioactive decay energy was measured in a scintilla- 
tion counter (Tri-CARB 2100 TR, Packard, Meriden, USA). 
The concentration of chl a+b in the extract was analysed 
using 96 % ethanol and used with chl 0+b concentrations 
from pigment analyses to express activities of Rubisco 
on a leaf surface area basis (jjimol C itT 2 s~ a ). Rubisco 
activation state (the percentage of Rubisco that was 
active in the leaves) was calculated as the ratio of 
initial activity to total activity. 



Chlorophyll fluorescence 

The maximum photochemical yield of photosystem II 
(PSII) in dark-incubated leaves, F v /F m , was measured in 
three dark-acclimated leaves from each genotype with 
a Portable Chlorophyll Fluorometer (PAM-2000, Walz 
Mess- und Regeltechnik, Effeltrich, Germany). Leaves 
were darkened with a leaf clamp for 15 min prior to 
measurements, and rapid light curves were measured 
using pre-installed software to estimate the relative effi- 
ciency of photon conversion (quantum yield) at low irra- 
diances (4>pam)> the light saturation point (I|<,pam, 
ixmol itT 2 s~ a ) and the maximum electron transport 
rate (ETR max , ixmol rrT 2 s~ a ). Differently from J max , also 
defined as electron transport rate and related to RuBP 
regeneration, ETR max refers to the electron transport in 
PSII activated by the light harvested by the chlorophylls. 

Statistics 

The software Statgraphics Centurion XVI (Statpoint Tech- 
nology Inc., VA, USA) was used to analyse the data. Data 
were tested for normal distribution and variance homo- 
geneity using Levene's test. The differences among the 
different P. australis haplotypes were tested by nested 
analysis of variance (ANOVA), with genotype nested in 
haplotype and using Type III sum of squares and the 
General Linear Model procedure of the Statgraphics soft- 
ware. Post hoc comparisons of means were performed 
using Tukey's honestly significant difference (HSD) pro- 
cedure at the 0.05 significance level. When necessary, 
data were log-transformed to approximate normality 
and secure variance homogeneity, but for clarity all 
data are presented as untransformed. A principal com- 
ponent analysis (PCA) was carried out to obtain a small 
number of linear combinations of parameters which 
vary significantly among haplotypes and genotypes, 
and account for most of the variability in the data. The 
components with eigenvalues >1 were extracted, and 
biplots of the principal components were obtained to 
illustrate differences among haplotypes and genotypes. 

Results 

Photosynthetic light responses 

The light-response curves of the four haplotypes differed 
significantly for all estimated parameters except the ap- 
parent quantum yield, (Table 2 and Fig. 1). The A max 
and I k varied significantly among genotypes within haplo- 
types. Haplotype M (including the EU-type and Greenyl- 
type) had consistently lower A max and I k than haplotype 
12 (the Land-type) and Ml (the Delta-type). The I c and 
Rdark of haplotype M were also lower than those of haplo- 
type 12 but did not differ from haplotype Ml. Haplotype Al 
(the Greeny3-type) did not differ in A max , I k , I c and R dark 
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Table 2. Photosynthetic characteristics of four haplotypes of Phragmites co-existing at the Gulf Coast of North America — light response 
parameters. Mean ( + 1 SE) light-saturated rate of photosynthesis (A max ), apparent quantum yield (<!>]), light compensation point (I c ), light 
saturation point (I k ), dark respiration (Rdark) and intrinsic water-use efficiency at maximum light intensity (WUEj) and the results (P values) of 
nested two-way ANOVA. Different superscript letters in the columns indicate significant differences [post hoc Tukey's HSD test, P < 0.05) among 
haplotypes. Analysis of variance results in bold indicate P values <0.05. 



Haplotype A max (|junol <l>i (mol C0 2 I c Ik Rdark (nmol WUEj (ixmol 

C0 2 m _2 s _1 ) (mol photons)" 1 ) (^mol m" 2 s" 1 ) (^mol m" 2 s" 1 ) C0 2 m" 2 s" 1 ) C0 2 mor 1 H 2 0) 



12 


36.1 + 2.3 b 


0.048 : 


h 0.002 


25.8 + 


1.7 b 


784 + 


42 b 


1.23 + 


0.09 b 


40.7 + 


2.3 


Ml 


34.3 + 2.0° 


0.052 : 


b 0.002 


19.6 + 


15 cb 


683 + 


36 b 


1.02 + 


0.08 ab 


51.6 + 


7.2 


M 


22.4 + 2.3 a 


0.049 : 


b 0.002 


15.1 + 


1.7 a 


475 + 


42° 


0.74 + 


0.09 a 


57.4 + 


6.0 


AI 


29.1 + 4.0 ab 


0.051 : 


b 0.003 


27.3 + 


3.0 ab 


602 + 


73 ab 


1.38 + 


0.16 ab 


41.4 + 


11.1 


Nested two-way ANOVA 


results (P value) 






















Genotype (haplotype) 


0.022 


0.283 




0.613 




0.043 




0.589 




0.002 




Haplotype 


0.030 


0.516 




0.022 




0.014 




0.036 




0.697 





30 




Figure 1. Average photosynthetic light-response curves of four haplotypes of Phragmites co-existing at the Gulf Coast of North America. The 
curves of haplotype 12 (Land-type), haplotype Ml (Delta -type), haplotype M (EU-type and Greenyl-type) and haplotype AI (Greeny3-type) are 
the average of three, four, three and one genotype, respectively (three light-response curves on three different leaves were analysed for each 
genotype). All the curves were conducted at 400 ppm C0 2 and 25 °C. 



either from haplotype M or 12 and Ml. The WUE, did not 
differ among the four haplotypes but differed among gen- 
otypes within haplotypes (Table 2). 

Photosynthetic C0 2 responses 

The parameters derived from the C0 2 -response curves 
varied among and within haplotypes (Table 3 and 
Fig. 2) except for the rate of mitochondrial respiration, 
R d , which did not differ among haplotypes. Significant 
differences were found for Vcmax and J max among haplo- 
types. The V cmax and J max were lowest in haplotype M 



and highest in haplotypes 12 and Ml. Values for haplo- 
type AI were intermediate and did not differ from 
those of haplotypes M and 12 and Ml. 

SLA and pigments 

Specific leaf area was higher in haplotype M than haplo- 
type 12. Both haplotypes M and 12 were not different 
from haplotypes Ml and AI (Table 4). Chlorophyll and 
total carotenoid contents did not vary among and 
within haplotypes. However, when pigment contents 
were expressed on an area basis (g m~ 2 ) they differed 
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Table 3. Photosynthetic characteristics of four haplotypes of Phragmites co-existing at the Gulf Coast of North America — C0 2 response 
parameters. Mean ( + 1 S.E) mitochondrial respiration (R d ), maximum carboxylation rate (V cmax ) and maximum electron transport rate (J max ) 
estimated from C0 2 -response curves of individual leaves and the results (P values) of nested two-way ANOVA. Different superscript letters in the 
columns indicate significant differences (post hocTukey's HSD test, P< 0.05) among haplotypes. Analysis of variance results in bold indicate 
P values <0.05. 



Haplotype 


R d ((jimol C0 2 m" 2 s" 1 ) 


Vcmax (fJ.mol m" 2 s" 1 ) 


Jmax (nmol m -2 s _1 ) 


12 


2.37 + 0.18 


118 + 6 b 


178 + 10 b 


Ml 


2.59 + 0.15 


123 + 6 b 


172 ± 9 b 


M 


2.27 + 0.18 


84 + 6 a 


116 + 10° 


AI 


1.97 + 0.31 


99 + ll ab 


137 + 17 ab 


Nested two-way ANOVA results (P value) 






Genotype (haplotype) 


0.006 


0.041 


0.014 


Haplotype 


0.411 


0.021 


0.019 
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Haplotype M1, Delta-type 




Haplotype Al, Greeny3-type 



/ 
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Intercellular CO2 concentration, Cj (umol mol 1 ) 

Figure 2. Average C0 2 -response curves of four haplotypes of Phragmites co-existing at the Gulf Coast of North America. The curves of haplo- 
type 12 (Land-type), haplotype Ml (Delta-type), haplotype M (EU-type and Greenyl-type) and haplotype AI (Greeny3-type) are the average of 
three, four, three and one genotype, respectively (three C0 2 -response curves on three different leaves were analysed for each genotype). As- 
similation values obtained with intercellular C0 2 concentration (CJ <250 (fjimol mol" 1 ) were fitted to the A c curve and assimilation values 
>250 (^imolmor 1 ) were fitted to the Aj curve (Dubois et al. 2007). All the curves were conducted at a light intensity of 
2000 |j.mol m" 2 s" 1 and 25 °C. 
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significantly among haplotypes, with lower contents in 
haplotype M than in haplotypes 12 and Ml. Haplotype 
AI did not differ from either haplotype M or 12 and Ml 
in pigment content (g rrT 2 ). 

Stomatal density and guard cell length 

Stomatal density and guard cell length on the adaxial 
and abaxial sides of the leaves did not differ among 
haplotypes (Table 5). However, there was significant vari- 
ation in adaxial stomatal density and abaxial guard cell 
length among genotypes within haplotypes. In general, 
leaves with longer stomata had lower stomatal densities. 

Rubisco activity 

The initial Rubisco activity, the total Rubisco activity and 
the Rubisco activation state did not differ among haplo- 
types. Only total Rubisco activity differed among geno- 
types within haplotypes (Table 6). The activity levels 
(24.1-32.0 fjimol C itT 2 s~ a for total Rubisco activity) 
(Table 6) were at the same level, or slightly lower, than 
the light-saturated gross rates of photosynthesis, A max 
(22.4-36.1 ixmol C0 2 nrT 2 s" 1 ) (Table 2). 

Chlorophyll fluorescence 

Chlorophyll fluorescence measurements revealed signifi- 
cant differences among genotypes within haplotypes for 
F v /F m and for I|<,pam from the rapid light-response curves 
(Table 7). The variability among genotypes within haplo- 
types was higher than the variability among haplotypes. 
Hence, there were no differences among haplotypes. 

Principal component analysis 

Of the significantly different parameters among haplo- 
types and genotypes, the PCA analysis extracted three 
components with eigenvalues >1. Together, they 
accounted for 88 % of the variability in the original data. 
The parameters contributing the most to PCI were A max , 
Vcmax, J max, 4, Rdark, total carotenoids and chl 0+b concen- 
trations expressed on a leaf-area basis (all positive load- 
ings) and SLA (negative loading). PC2 had high positive 
loadings from WUE, and R d , and negative loadings for 
adaxial stomatal density, total Rubisco and I|<,pam- PC3 
had high positive loadings of abaxial stomata guard cell 
length, I c and F v /F m . The biplots of PCI vs. PC2 (Fig. 3A) 
and PCI vs. PC3 (Fig. 3B) separated the genotypes of the 
four haplotypes. Haplotypes 12 and Ml showed higher 
PCI values than haplotype M (Fig. 3A). Haplotype 12 was 
separated from haplotype Ml along the PC3 axis 
(Fig. 3B). Haplotype AI was located in the lower part of 
the biplot (Fig. 3A), reflecting its relatively low PC2 
scores. Also, ROMS7 (haplotype Ml, Delta-type) was dis- 
tinct from the cluster of haplotype Mi's genotypes along 
the PC2 axis (Fig. 3A). ROMS7, unlike the other genotypes 
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Table 5. Characteristics of four haplotypes of Phragmites co-existing at the Gulf Coast of North America — stomatal density and guard cell 
length. Mean ( + 1 SE) stomatal density and guard cell length on the adaxial and abaxial leaf surface and results (P values) of nested two-way 
ANOVA. There are no significant differences among haplotypes (post hoc Tukey's HSD test, P < 0.05). Analysis of variance results in bold indicate 
P values <0.05. 



Haplotype Adaxial density (mm 2 ) Abaxial density (mm 2 ) Adaxial length (|jim) Abaxial length (i-im) 



12 


631 + 42 


681 + 


50 


22.1 + 


0.7 


22.7 + 


0.7 


Ml 


666 + 37 


644 + 


43 


20.1 + 


0.6 


19.9 + 


0.6 


M 


664 + 42 


712 + 


50 


20.2 + 


0.7 


19.6 + 


0.7 


AI 


793 + 74 


822 + 


86 


19.8 + 


: 1.2 


20.2 + 


1.3 


Nested two-way ANOVA 


results (P value) 














Genotype (haplotype) 


0.023 


0.158 




0.068 




0.012 




Haplotype 


0.438 


0.419 




0.419 




0.097 





Table 6. Characteristics of four haplotypes of Phragmites co-existing 
at the Gulf Coast of North America — Rubisco activity. Mean ( + 1 SE) 
initial Rubsico activity (non-activated Rubisco), total Rubisco activity 
and Rubisco activation state (% of active Rubisco) and results 
(P-values) of nested two-way ANOVA. There are no significant 
differences among haplotypes (post hoc Tukey's HSD test, P < 0.05). 
Analysis of variance results in bold indicate P values <0.05. 



Haplotype 


Non-activated 


Total Rubisco 


Rubisco 




Rubisco (ixmol 


((imol 




activation 




Cm" 2 s"') 


Cm" 2 


s" 1 ) 


state (%) 


12 


14.3 + 1.8 


26.8 + 


2.8 


53.9 + 2.1 


Ml 


15.9 + 1.5 


24.1 + 


2.4 


64.2 + 1.9 


M 


16.2 + 1.8 


27.4 + 


2.8 


59.1 + 2.1 


AI 


21.2 + 3.2 


32.0 + 


4.9 


66.0 + 3.7 


Nested two-way ANOVA results (P value) 






Genotype 


0.132 


0.027 




0.852 


(haplotype) 










Haplotype 


0.454 


0.598 




0.055 



of haplotype Ml, shares alleles with the EU-type geno- 
types of haplotype M (allele 202 at locus PaGT22; Lamber- 
tini et a/. 2012a). Within haplotype M, the Greenyl-type did 
not appear different from the EU-type and was distinct 
from the Greeny3-type with haplotype AI along the PC2 
axis. 

Discussion 

We observed substantial differences among haplotypes 
and genotypes within haplotypes in the Gulf Coast of 
North America. Haplotypes 12 (Land-type) and Ml 
(Delta-type) had higher A max , V cmax , J max , I k , leaf carot- 
enoids and chl a+b concentrations expressed on a 



leaf-area basis than haplotype M (EU and Greenyl- 
types) (Tables 2-4). Haplotype AI had intermediate 
values for these parameters and was not different 
from either haplotype M or 12 and Ml. 

Genotypes varied significantly in WUE,, R d > adaxial sto- 
matal density, abaxial stomatal guard cell length, total 
Rubisco activity, F v /F m and Ir.pam- However, as genotypes 
clustered according to haplotype in the PCAwhen signifi- 
cant variation among genotypes was considered, such 
differences appear to be within the normal variation 
range within haplotypes. Only one genotype (ROMS7), 
probably a hybrid between haplotype Ml and haplotype 
M (Lambertini et al. 2012a), deviated within the cluster of 
haplotype Mi's genotypes and showed higher WUE, and 
mitochondrial respiration (R d ), and lower I|<,pam and 
adaxial stomatal density, than its close relatives within 
haplotype Ml and haplotype M. The geographic range 
in which gene flow between these two haplotypes 
occurred (Gulf Coast vs. Mediterranean region) is, 
however, still to be defined (Lambertini et al. 2012a, b), 
but the genetic and ecophysiological variation within 
haplotype Ml with the Delta-type phenotype is present 
in the Gulf Coast and deserves to be studied further in re- 
lation to the phylogenetic relationships of this group. 
Within haplotype M, the Greenyl-type genotype (OCT1) 
appears to be more similar in photosynthetic traits to 
the EU-type genotypes sharing haplotype M than to 
the genotype with haplotype AI sharing the Greeny 
phenotype. 

As the differences among genotypes have a genetic 
basis — and are thus inheritable — the different haplotypes 
can be regarded as ecotypes. The photosynthetic charac- 
teristics of the haplotypes, their recent introduction to the 
Gulf Coast (Lambertini et al. 2012a), as well as the fact 
that they grow in direct contact (adjacent stands) in the 
Mississippi River Delta, suggest that the adaptations 
were acquired in the native, rather than in the introduced 
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Table 7. Characteristics of four haplotypes of Phragmites co-existing at the Gulf Coast of North America — chlorophyll fluorescence. Mean ( + 1 
SE) potential quantum yield of PSII (F v /F m ), quantum yield (c&pam). maximum electron transport rate (ETR max ) and light saturation point (Ik.pam) 
of the rapid light-response curves from chlorophyll fluorescence and the results (P values) of nested two-way ANOVA. There are no significant 
differences among haplotypes (post hoc Tukey's HSD test, P< 0.05). Analysis of variance results in bold indicate P values <0.05. 



Haplotype 


FJF m 




'I'PAM 




ETR max (.xmol m" 2 s" 1 ) 


Vam (t*m°l m~ 2 s" 1 ) 


12 


0.796 + 


0.005 


0.319 


+ 0.006 


449 + 74 


142 + 30 


Ml 


0.788 + 


0.005 


0.311 


+ 0.005 


599 + 64 


164 + 26 


M 


0.809 + 


0.005 


0.318 


± 0.006 


294 + 74 


110 + 30 


AI 


0.798 + 


0.009 


0.285 


+ 0.010 


798 + 129 


227 + 51 


Nested two-way ANOVA results (P value) 












Genotype (haplotype) 


0.048 




0.131 




0.191 


0.026 


Haplotype 


0.136 




0.079 




0.055 


0.267 



A 4.7 
2.7 

a 07 

a. 

-1.3 
-3.3 

B 3.3 
2.3 
1.3 

u 

Q_ 

0.3 
-0.7 
-1.7 



□ ROMS? 



ROM16 

AtWttJr Scrruu guard cell length- 



OCTKGfeaflyi) 




TotJRubiseG y DPalOTgcUStCifMnyi ' 



-6.3 



-4.3 



-2.3 -0.3 
PCI (51.7%) 



1.7 



3.7 



Haplotype 12, Land- type 
Haplotype Ml, Delta-type 
. Haplotype M, EU- and Greenyl -types 

Haplotype Al, Greeny3-type 



MUM 
□ 

ADariit stomm guard cell length 




PCI 



Figure 3. Biplots from a PCA based on all significantly different ecophysiological parameters among haplotypes and genotypes of Phragmites 
co-existing at the Gulf Coast of North America. PCI, PC2 and PC3 explain 51.7, 25.7 and 10.6 % of the variation, respectively. 



range. Nguyen et al. (unpubl. data) also observed the 
same ecophysiological differences among genotypes of 
the same haplotypes sampled in the native range. 

The African and Mediterranean haplotypes (12 and 
Ml) are adapted to warmer climatic conditions than 



haplotype M, which has photosynthetic characteristics 
more similar to European P. australis in the temperate 
region (Lessmann et al. 2001; Hansen et al. 2007). 
Salvucci and Crafts-Brandner (2004) found similar 
results for plants from the desert, like Larrea tridentata of 
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the southwest USA, vs. Deschampsia antarctica, endemic to 
maritime Antarctica. The other European haplotype AI had 
intermediate photosynthetic characteristics between the 
Eurasian haplotype M and the African and Mediterranean 
haplotypes 12 and Ml. Phytogenetically, haplotype AI is 
closely related to a P. australis population in South Africa, 
Namibia and Botswana, a population that is introduced 
to Europe (Lambertini et al. 2012a, fa). 

In the present study, the A max of haplotype M 
(22.4 |jimol nrT 2 s~ a ) was within the range reported for 
P. australis grown in Europe (11.9-23.3 jjimol rrT 2 s~ a ) 
(Hansen et al. 2007) and in the introduced range in 
North America (19.2-22.8 jjimol) (Mozdzer and Zieman 
2010). The A max of haplotype AI (Greeny3-type) 
(29.1 + 4.0 |j,mol rrT 2 s~ a ) was beyond the upper range 
of A max values in these European populations (Hansen 
et al. 2007). The V cmax and J max of the four haplotypes 
were consistent with A max , which were higher in haplo- 
types 12 and Ml than in haplotype M, but not different 
from haplotype AI. 

In our study, the A max of haplotype M was lower than 
those of haplotypes 12 and Ml (Land- and Delta-types: 
36.1 and 34.3 |j,mol rrT 2 s~ a ). In dark conditions, haplo- 
type 12 respired more than haplotype M, as expected by 
the model of Prioul and Chartier (1977) for plants with 
higher A max . However, R d ari< of haplotype Ml did not 
differ from that of the other haplotypes and had 
higher A max , I c and I k , i.e. typical traits of a sun-adapted 
plant. This light response is expected for plants from 
warm and arid regions (Sage and Monson 1999). 

The leaves of haplotype 12 were thicker (lower SLA) 
than those of haplotype M, and leaves of haplotypes 12 
and Ml also contained higher contents of photosynthetic 
pigments than haplotype M, when expressed on a leaf 
surface area basis, probably due to more chloroplasts of 
smaller size in the thick leaves (Demmig-Adams and 
Adams 1992; Lambers et al. 2008). This is also consistent 
with sun plants (Larcher 2003). A higher number of chlor- 
oplasts in the leaves increases the number of Calvin 
cycles and thus the photosynthetic rates. 

The studied haplotypes grow in direct contact in the 
Mississippi River Delta, and hence compete, but their 
broad distribution in the delta may also be determined 
by factors such as hydrological regime and salinity 
(Vasquez et al. 2005). The introduced Eurasian 
P. australis that has invaded widely in North America, 
and which constitutes the source population of the 
EU-type phenotype of haplotype M, has been studied in- 
tensively. This introduced haplotype M has been spread- 
ing and outcompeting the native P. australis subsp. 
americanus (Saltonstall et al. 2004) by means of faster 
growth, greater investment in sexual reproduction and 
more efficient seed dispersal (Howard et al. 2008; 



Belzile et al. 2010; McCormick et al. 2010, Kettenring 
et al. 2011). Furthermore, this introduced haplotype M 
also tolerates higher salinity levels (Vasquez et al. 
2005), and has a greater ventilation efficiency for 
rhizome aeration than the native P. australis ssp. 
americanus (Tulbure et al. 2012). 

Our study suggests that haplotypes 12 and Ml are 
superior to haplotype M in terms of photosynthetic C0 2 
assimilation. The unique genotype of haplotype AI 
analysed in this study lies between haplotypes 12, Ml 
and M. This is consistent with the current distribution 
of the P. australis haplotypes in the Mississippi Delta, 
with haplotype 12 (Land-type) dominating upland and 
haplotype Ml (Delta-type) dominating in the delta 
marshes. The different distributions of these two domin- 
ating haplotypes may be due to different ecophysiologi- 
cal characteristics not investigated in this study, e.g. 
differences in salinity tolerance and/or differences in 
the ability to aerate rhizomes or roots. Environmental 
differences between the delta and the uplands may 
also be important. Sediment accretion, variable water 
levels and, not least, tropical hurricanes provide 
windows of opportunity for the establishment of 
P. australis seedlings more frequently in the delta than 
in the progressively consolidated lands upstream of the 
delta (Cretini et al. 2012). Haplotypes M (both EU- and 
Greenyl-types) and AI (Greeny3-type) are sympatric in 
the delta with haplotype Ml, but appear less proficient 
in photosynthetic C0 2 assimilation than haplotype Ml. 
However, these phenotypes seem to be spreading at 
the expense of the Land- and Delta-types (Hauber 
et al. 2011). This discrepancy might be explained by 
several factors. First, photosynthetic performance is 
only one of many factors determining the competitive 
ability of the P. australis haplotypes in the delta. Physical 
disturbances, water regime, salinity (Vasquez et al. 2005) 
and biotic factors such as resistance to grazing and 
various pests are likely to be very important (Blossey 
and Noetzold 1995). Secondly, our study was performed 
on plants in a controlled steady-state environment. In 
nature, the prevailing environmental conditions are 
changing on different temporal scales, and it might be 
extreme conditions, rather than the average conditions, 
that determine the competitive ability, and hence the 
fitness, of a plant (Wilcove et al. 1998). The origins of 
the four haplotypes and adaptations derived from their 
native climatic zones nevertheless seem to explain the 
differences in their photosynthetic ability. 

Conclusions 

Our study has documented that four P. australis haplo- 
types co-existing in the Mississippi River Delta differ in 
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important photosynthetic traits, and we suggest that 
these differences in photosynthetic characteristics are 
related to the climatic conditions in the native range of 
the haplotypes. It remains, however, to be understood 
how these differences, ecologically definable as 'photo- 
synthetic ecotypes', relate to plant competitiveness 
and fitness in the Gulf Coast environment. Another 
aspect that deserves further attention is intra-haplotypic 
variation. Our study shows that hybrids might differ in 
ecophysiological traits— and hence in adaptive strat- 
egies — from both their parents, and that genotypes 
showing the same phenotype can have different origins, 
as well as different ecophysiological characteristics. 
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